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ISOPHOT 95 fim observations in the Lockman Hole 

The catalogue and an assessment of the source counts* 



o 
o 

(N 

Sh 

Oh 
< 

in 



> 
o 

m 

o 

o 

Or 

o ■ 
& . 



G. Rodighiero 1 and A. Franceschini 1 

Dipartimento di Astronomia, Universita di Padova, Vicolo dell'Osservatorio 2, 1-35122 Padova, Italy e-mail: 
rodighiero@pd . astro . it 

Received ; accepted 

Abstract. We report results from a new analysis of a deep 95 pm imaging survey with ISOPHOT on board the Infrared Space 
Observatory, over a ~1 deg 2 area within the Lockman Hole, which extends the statistics of our previous study (Rodighiero 
et al. 2003). Within the survey area we detect sixty-four sources with S/N> 3 (roughly corresponding to a flux limit of 16 
mjy). Extensive simulations indicate that the sample is almost complete at Sgsim > 100 mjy, while the incompleteness can be 
quantified down to ~30 mjy. The 95 pm galaxy counts reveal a steep slope at S 95,,,,, < 100 mly (a = 1.6), in excess of that 
expected for a non-evolving source population. In agreement with counts data from ISO at 15 and 175 /jm, this favours a model 
where the IR populations evolve both in number and luminosity densities. We finally comment on some differences found with 
other ISO results in this area. 
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1. Introduction 

The infrared sky observed from space is providing tight con- 
straints on the evolution of cosmological source populations. 
The IRAS extragalactic number counts (Oliver et al. 1992, 
Bertin et al. 1997) showed some marginally significant excess 
of galaxies compared to non-evolutionary predictions. More re- 
cently, ISOPHOT (Lemke et al. 1996) found some evidence of 
evolution for the IR population at 170/im(e.g. Pugetetal. 1999, 
Dole et al. 2001), but source confusion at these wavelengths 
limits the observations to moderately faint fluxes (5 175^,,, > 135 
mjy). At shorter wavelengths, ISOCAM (Cesarsky et al. 1996) 
detected a number of faint sources consistent with strong evo- 
lution in the mid-IR (e.g. Elbaz et al. 1999), a factor ~10 in 
excess with respect to non-evolution predictions. 

The ISO deep counts, together with those in the submil- 
limeter (e.g. Scott et al. 2002) are used as constraints on mod- 
els of galaxy formation and evolution (e.g. Lacey et al. 1994, 
Franceschini et al. 2001, Lagache et al. 2003). These studies 
are also needed to estimate the confusion noise for future IR 
and submillimeter telescopes (Xu et al. 2001, Lagache et al. 
2003). 

In this paper we present a new analysis of ISOPHOT 95 
pm data in the direction of the Lockman Hole, extending our 
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previous study (Rodighiero et al. 2003, hereafter Paper I) and 
significantly improving the source statistics. 

The Spitzer Space Observatory (Fazio et al., 1999) will 
soon observe the same region of the Lockman Hole in com- 
plementary wavebands at 24, 70 and 160 /mi with MIPS and 
in the near-IR with IRAC. This will provide additional exten- 
sive information on the Spectral Energy Distributions (SEDs) 
of ISO sources. 

The paper is organized as follows: in Section 2 we intro- 
duce the dataset and in Section 3 we summarize features of the 
adopted reduction procedure. The infrared maps and the new 
catalogue are presented in Section 4. The source counts and a 
discussion of the effects of confusion and cosmic variance are 
reported in Section 5. 

2. The sample: LHEX + LHNW 

The Lockman Hole (Lockman et al., 1986) is the sky region 
with the lowest HI column density, and is thus best suited for 
the detection of faint extragalactic infrared sources. Moreover, 
available multiwavelength observations in this area (Hasinger 
et al. 2001, Fadda et al. 2004 in prep., and Rodighiero et al., 
2004 in prep., Rodighiero et al., 2003, Kawara et al., 2004, 
Scott et al., 2002, De Ruiter et al., 1997, and Ciliegi et al. 2003), 
allow to study the spectral energy distributions of the detected 
sources and to address their physical nature. 

ISOPHOT, on board the ISO satellite, observed two differ- 
ent regions in the Lockman Hole. Each one of the two fields, 
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3. Data reduction 
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The LHNW data set has been reduced with the procedure de- 
scribed in Paper I and successfully applied to the LHEX field. 
The procedure is based on a parametric algorithm that fits the 
signal time history of each detector pixel. It extracts the back- 
ground level, in order to identify the singularities induced by 
cosmic-rays impacts and by transient effects in the detectors, 
and to identify and extract real sky sources. 

The source extraction process, the completeness, the pho- 
tometric and astrometric accuracies of the final catalogues have 
been tested by us with extensive sets of simulations, by insert- 
ing into the real image sources with known flux and position, 
and accounting for all the details of the reduction procedure. 
We defer the interested reader to Paper I for details. The flux 
calibration was also verified by reducing with the same tech- 
nique ISOPHOT 95 fim observations of some selected standars, 
such as calibrating stars and IRAS sources, in order to cover the 
flux range spanned by our Lockman source catalogue. 
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Fig. 1. The 95 fim map of the Lockman Hole LHEX (upper 
panel) and LHNW (lower panel) areas. The circles indicate our 
detected sources with signal-to-noise ratio > 4. Each map is 
~ 44' X 44'. 



called LHEX and LHNW, covers an area of ~ 44' x 44' and has 
been surveyed at two far-infrared wavebands with the CI 00 and 
C200 detector (respectively at 95 and 175 //m), in the P22 sur- 
vey raster mode. The ISOPHOT survey in the LHEX field has 
been reported in Paper I, while we extend here our analysis to 
the LHNW, which consists of a mosaic of four rasters, each one 
covering an area of ~ 22' x 22'. The observational parameters 
are reported in Paper I. 



4. The 95 /mi source catalogue 

In Paper I, we detected 36 sources with S/N> 3 (correspond- 
ing to a flux of 16 mJy) within the LHEX area, making up a 
complete flux-limited sample for Sg^ m > 100 mJy. Reliable 
sources were detected, with decreasing but well-controlled 
completeness, down to Sgs^m — 20 mJy. In the twin LHNW 
area of ~0.5 deg 2 , we detect 28 sources above a 3-cr threshold, 
reaching the same flux level as in the previous sample. 

In Figure [Qwe show the final maps obtained in the LHEX 
and in the LHNW fields. The open circles indicate sources de- 
tected with a signal-to-noise ratio greater than 4. In Table ^ 
we report the new list of sources detected in the LHNW re- 
gion. The entries are: IAU-conformal names, sky coordinates 
(right ascension and declination at Equinox J2000), the detec- 
tion significance (signal-to-noise ratio), the 95 fim total fluxes 
(in mJy) and their uncertainties. As discussed in Paper I, all 
sources have been extracted from the maps and confirmed by 
visual inspection of the pixel history. This approach produces 
an highly reliable source catalogue. 

Our source extraction tool and simulation procedures can 
be optimally used to recover the total fluxes for slightly blended 
sources. An example is reported in Figure [2] where for a pair 
of ISOPHOT complex sources (sources 8 and 17 in LHEX, 
see catalogue in Paper I) we show the optical R band image 
(Fadda et al. 2004, in prep.) with overlayed the 95 fim con- 
tours. The ISOPHOT positions are marked with open squares, 
the triangles indicate the ISOCAM 15 fim detections (Fadda et 
al., 2004). The distance between the two far-IR sources, ~ 55 
arcseconds, typical for the few blended objects in the LHEX 
and LHNW maps, slightly exceeds the ISOPHOT CI 00 beam 
(~ 45 arcsec). 

The final combined LHEX+LHNW catalogue contains a 
total of 64 sources, detected at 95 fim over an area of about 1 
square degree. 
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Table 1. LHNW source catalogue. 



ID 


RA 


DEC 


S/N 


Flux 




(J2000) 


(J2000) 




TmJvl 


LHJ103521 +580034 


10:35:21.2 


+58:00:34 


33 


235 ± 24 


LHJ103606+574715 


10:36:06.7 


+57:47:15 


22 


163 ± 17 


LHJ103604+574815 


10:36:04.4 


+57:48:15 


21 


181 ± 20 


LHJ103653+575444 


10:36:53.8 


+57:54:44 


16 


201 ± 24 


LHJ 103618+573904 


10:36:18.2 


+57:39:04 


16 


227 ± 26 


LHJ103515+573330 


10:35:15.5 


+57:33:30 


16 


136 ± 16 


LHJ 103523+575656 


10:35:23.3 


+57:56:56 


15 


114 ± 14 


LHJ103318+574925 


10:33:18.6 


+57:49:25 


15 


159 ± 19 


LHJ103236+588042 


10:32:36.2 


+58:80:42 


12 


109 ± 14 


LHJ 103537+573537 


10:35:37.4 


+57:35:37 


12 


114 ± 13 


LHJ103530+573139 


10:35:30.8 


+57:31:39 


8 


73 ± 12 


LHJ 103409+5727 1 5 


10:34:09.1 


+57:27:15 


8 


63 ± 10 


LHJ103407+572613 


10:34:07.5 


+57:26:13 


8 


69 ± 10 


LHJ103341+582031 


10:33:41.2 


+58:20:31 


8 


77 ± 12 


LHJ103202+588009 


10:32:02.4 


+58:80:09 


8 


97 ± 15 


LHJ103610+574330 


10:36:10.6 


+57:43:30 


7 


85 ± 14 


LHJ103336+582000 


10:33:36.7 


+58:20:00 


6 


34 ± 7 


LHJ 103308+588030 


10:33:08.8 


+58:80:30 


5 


80 ± 16 


LHJ 103445+574733 


10:34:45.1 


+57:47:33 


5 


48 ± 10 


LHJ103249+573719 


10:32:49.7 


+57:37:19 


5 


56 ± 1 1 


LHJ 103352+584030 


10:33:52.1 


+58:40:30 


4 


55 ± 14 


LHJ 103304+573637 


10:33:04.3 


+57:36:37 


4 


55 ± 14 


LHJ 103550+573005 


10:35:50.4 


+57:30:05 


4 


49 ± 11 


LHJ103159+587026 


10:31:59.8 


+58:70:26 


4 


100 ±24 


LHJ103524+575749 


10:35:24.5 


+57:57:49 


3 


35 ± 7 


LHJ103316+573136 


10:33:16.1 


+57:31:36 


3 


30 ± 10 


LHJ103327+574539 


10:33:27.6 


+57:45:39 


3 


42 ± 12 


LHJ1033 13+589010 


10:33:13.2 


+58:90:10 


3 


21 ± 7 




Fig. 2. Example of a pair of ISOPHOT complex sources. The 
Figure shows the optical R band image (Fadda et al. 2004, in 
prep.) with overlayed the 95 /im contours. The ISOPHOT po- 
sitions are marked with open squares, the triangles indicate the 
ISOCAM 15 jum detections (Fadda et al., 2004). 



5. Source counts from the Lockman Hole 95 //m 
survey 

We have computed the 95 fim source counts down to a flux level 
of 30 mJy, following the same procedure described in Paper I. 
The counts have been corrected for the incompleteness factors 
derived through simulations. According to these, the complete- 
ness decreases from ~ 80% around 100 mJy, to ~ 30% around 
40 mJy, see Table 1 in Paper I). Poissonian noise has been con- 
sidered. 

The new source counts confirm our previous finding for 
the LHEX field. The integral counts are plotted in Figure 
[3] as starred symbols. Our results are compared with those 
from other surveys: the preliminary analysis of the ISOPHOT 
ELAIS survey (Efsthatiou et al. 2000, open circles); the counts 
derived from the IRAS lOOyum survey (open squares); those 




.000 0.100 0.010 

Flux Density [Jy] 



Fig. 3. Integral counts at 95 fan. Our estimated values (starred 
filled symbols) are compared with data from a variety of other 
surveys (see text for references and details). The counts are 
expected to be dominated by extragalactic sources, for which 
some model predictions are also reported. 



published by Linden- Vornle et al. (2000, open triangles), Juvela 
et al. (2000, '+' symbol) and recently by Kawara et al. (2004, 
filled triangles) which are based on the same dataset (see dis- 
cussion below). A comparison is made also with model counts 
by Franceschini et al. (2001, solid line), Lagache et al. (2002, 
dot-dashed line), Xu et al. (2001, dashed line) and Rowan- 
Robinson et al. (2001, dotted line). The slope of the counts 
turns out to be a ~ 1.6. 

The good agreement with the ELAIS source counts (which 
are based on a larger area) suggests that these results are not 
much affected by spatial biases. In the fainter flux region be- 
low 100 mJy, the source counts from different surveys appear 
instead more scattered. 

In Figure we report the differential 95 /mi counts dN/dS 
normalized to the Euclidean law /V oc S~ 2 5 (the data are re- 
ported in Table |3. For this we have used only sources with 
a signal-to-noise ratio greater than 4. The results are com- 
pared with those from Kawara et al. (2004, filled triangles) 
and from the preliminary analysis of the ISOPHOT ELAIS sur- 
vey (Efsthatiou et al. 2000, open circles). A comparison is per- 
formed with modellistic differential counts from Franceschini 
et al. (2001), Xu et al. (2001) and Lagache et al. (2002). 

Our results differ from those of Kawara et al. (2004), who 
analysed the same dataset used here (i.e. the combined LHEX 
and LHNW ISOPHOT 95 yum data), but with an independent 
reduction based on the standard PHOT Interactively Analysis 
(PIA, Gabriel et al. 1997). Their final catalogue includes quite 
more sources than ours (223 versus 64), at similar flux levels 
and S/N ratios. The two number counts turn out to be in good 
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Table 2. Integral and differential source counts at 95 fan. 



^Lockman (Kowara '04) 
□ IRAS 
Q ELAIS 

^Lockman (this work) 




_ Franceschini 
H Lagache 



10.00 



S ra ,„ [Jy] 



Fig. 4. Differential 95 yum counts dN/dS normalized to the 
Euclidean law {N cc S~ 2 5 ). Symbols are the same as in Fig. 
13 Our results are compared with those from the preliminary 
analysis of the ISOPHOT ELAIS survey (Efsthatiou et al. 
2000, open circles) and with those from Kawara et al. (2004, 
filled triangles). We also show as a long-dashed line the ex- 
pected contribution of non-evolving spirals as in the model of 
Franceschini et al. (2001). 



agreement only above 200 mjy. Although unable to understand 
the origin of this difference, we note that different flux calibra- 
tions were adopted by the two analyses, Kawara et al. using a 
single reference, IRAS F10507+5723. In any case, our careful 
detector time-sequence analysis and the simulations performed 
to check each sources should in principle allow us a robust re- 
jection of spurious detections, while at the same time keeping 
under control the survey incompleteness. 

With this new extended sample we confirm our conclusions 
in Paper I about the existence of an evolving population of 
IR galaxies, consistent with independent results based on deep 
ISOCAM mid-IR and ISOPHOT 175 /jm counts. Our analysis, 
however, indicates that the evolution cannot be too strong since 
the level of the counts at the faintest fluxes is not very high. The 
counts predictions reported in the figures typically assume that 
the whole local galaxy population evolves back in cosmic time 
in luminosity and/or number density. With such an assumption, 
the model counts tend to more or less exceed our data at the 
faint fluxes. In the evolutionary scheme by Franceschini et al. 
(continuous lines) the contribution of a non-evolving normal 
galaxy population dominating the bright counts adds to that of 
a fast evolving starburst-galaxy component rising up at faint 
fluxes, then producing relatively flat overall counts. The present 
results at 95 yum seem to confirm this model. 

An apparently non-random spatial distribution of the 95 
/im sources was already indicated in Paper I. A comparison 
of the two maps in LHEX and LHNW (Fig. ^ shows even 
more clearly the effects of field-to-field variation (28 versus 36 
sources detected in the to areas). The source spatial distribu- 
tion in LHNW seems particularly inhomogeneous, with source 
clustering close to two map borders. This strong clustering of 
the far-IR sources (whose physical origin may be related with 



Integral 




Differential 


S dN(>S) 


flux bin 


bin center 


dN IdS x S ' 1 * 


mjy deg" 2 


mjy 


mJy 


deg~ 2 mJy L5 


30 70.0 


30-100 


54 


2.63e4± 1.06e4 


60 45.7 


60-150 


94 


3.42e4± 1.1 2e4 


100 18.8 


100-300 


173 


3.79e4± 1.25e4 


150 10.3 








300 1.8 









the effects of galaxy interactions and mergers as a trigger of the 
starburst-induced far-IR emission) makes the confusion prob- 
lem very serious. Four sources in LHNW and five in LHEX 
appear with close companions (see one example in Fig|2jl, em- 
phasizing that source blending and confusion happens even at 
these moderately faint flux limits. Although progress is being 
achieved with Spitzer, definitely larger space telescopes (e.g. 
the Herschel Space Observatory, due to launch in 2007) will be 
required for proper sampling of this difficult waveband domain. 
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